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[Continued from page 48.] 

This surface (43) might be called the probability surface, just as (1) is the 
probability curve. Its form is such that if it is intersected by a plane par- 
allel to the plane of XY, the lines of intersection will be ellipses whose 
centre and axes, projected on the plane of XY, coincide with the origin and 
axes of X and Y, and the axes of the ellipse in the x and y directions are 
inversely proportional to h^ and h 2 , and therefore directly proportional to 
j-j and r 2 . This may be shown by making z a constant in (43). The pro- 
jected ellipses are the loci of points of equal probability. 

From this symmetry of the surface with respect to the axes of X and Y, 
it follows that these are the free axes of the system of expanded coefficients 
z, regarded as masses in the plane of XY. Thus the given probabilities 
),, and the resultant probabilities z, have both the same free axes. This also 
follows from the permanence of direction of the axes which give maximum 
and minimum values to the radii of gyration in successive powers, as al- 
ready noticed. The value of 2 is a maximum at the origin, and this being 
the place of the centre of gravity of the whole system of masses *, it appears 
that the point of most probable error in the result of a very large number 
of observations, will be located at the centre of gravity of all the points 
which can possibly result under the given system of elementary errors, each 
possible point of error being taken with a weight proportional to the proba- 
bility that such error will occur. This also follows from what we have 
already proved in connection with formula (3). In order that the most 
probable result should coincide with the true point or place of zero error, it 
would be necessary and sufficient that the probabilities or weights X of the 
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elementary points of error should be such that their centre of gravity falls 
at the zero point. 

Suppose the plane of ZF to be divided up into elementary rectangles 
dxdy by lines drawn parallel to the axes ; then (43) gives the probability z 
that an error which occurs will fall on the rectangle whose centre is at the 
point x = idx, y=jdy. It gives an approximate result when instead of dx 
and dy we use small finite distances Ax and Ay. The values of h 1 and h 2 are 
computed by (37) and (38) if the elementary errors and their probabilities X 
are known. But if these are unknown, as generally happens in practical 
applications, we must compute from a large number of observed errors, the 
quadratic mean errors in the x and y directions, answering to r 1 |/fc and 
r 2 j/& in (38), and thus obtain h 1 and /t 2 . For example, if a large number 
of shots have been fired and have struck a target, and we wish to infer, from 
the distribution of the shot-marks, what is the law of probable distribution 
of other shots to be made under like conditions, our investigation leads to 
the following process. Regarding each given shot-mark as the mass of a 
material point, all such masses being equal, find the centre of gravity of the 
system, and then by (31) find its free axes. The distance of the centre of 
gravity from any assumed axis of reference is the arithmetical mean of the 
distances of all the shot- marks from that axis, distances on opposite sides of 
it having contrary signs. The free axes will be at right angles to each oth- 
er at the centre of gravity, and we take either one of them as the X and the 
other as the Y axis. Find the quadratic mean errors in the x and y direc- 
tions, that is, the square root of the mean of the squares of the distances of 
the shot-marks from the Y and X axis respectively. Denoting these by 
e 1 dx and e$dy, we shall have as in (38), 

{\dxf = l+2ej, {h 2 dyf = l-2ei, (44) 

and (43) becomes 

2^e 1 e 2 * v ' 

For convenience, dx and dy may be taken equal to each other and to any 
small unit of measure, for instance one inch. Then (45) gives the proba- 
bility z that any particular future shot will strike that square inch of the 
target whose coordinates are i and j inches. 

All this is on the supposition that we know nothing about the probable 
distribution of the errors except what is to be inferred from the actual dis- 
tribution of the shot-marks already made. 

The principle that the axes of X and Y should be taken to coincide with 
the free axes of the given shot-marks, is one which has not, so far as I know, 
been brought out by any previous writer. The usual course has been, to 
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draw the axes horizontally and vertically through the centre of gravity of 
the marks. This, however, will in general be correct only when there is 
reason to believe that the elementary errors, whose probabilities are 1, are dis- 
tributed so as to make the hoinzontal and vertical axes free axes, which will 
be the case for instance when these errors and their probabilities fall symmet- 
rically on either side of the horizontal and vertical axes, making them axes 
of symmetry in the system of masses L Granting that the free axes are 
horizontal and vertical, formula (43) agrees with the one in general use. It 
can be written 

z _ /hjdx e _i lW \ (fhdy- hV \ ^ 

the second member being the product of two factors of the same form as in 
the probability curve (1). This property has enabled writers to give a de- 
monstration of (43), or what has seemed to be one, in a very simple manner. 
(See for instance Didion's Calcul des Probabilite s applique au Tir des Pro- 
jectiles, Paris, 1858, p. 42.) If we suppose the plane of XY to be divided 
by lines drawn parallel to the Y axis, into bands of the width dx, and as- 
sume that the probability of deviation from the Y axis is the same at all 
parts of that axis, or in other words, is independent of y, then the proba- 
bility that an error which occurs will fall within the band whose distance 
from the Y axis is x, is represented by the first of the two factors in (46). 
Likewise if lines are drawn parallel to the X axis, forming bands of the 
width dy, the probability that the error will fall within the band whose 
distance from the X axis is y, is represented by the second factor in (46). 
But the probability that two independent events will both occur, is the pro- 
duct of their separate probabilities, so that (46) gives the probability z that 
the error will fall within both bands, that is, at their intersection, and with- 
in the rectangle dxdy whose coordinates are x and y. This demonstration, 
however, has been objected to, on the ground that the independence of the 
x and y deviations, if it exists, is not so evident as to be axiomatic. It re- 
ally seems to be a thing to be proved. Various authors who have discus'd 
the point are named in Merriman's List of Writings relatiiig to the Method 
of Least Squares. (Trans. Conn. Acad., Vol. IV, Part 1.) Our present 
investigation shows that the x and y deviations may be regarded as inde- 
pendent of each other, provided that the axes of X and Y are taken to 
coincide with the free axes of the probabilities k of the elementary errors. 
When these errors are unknown, we approximate to the desired axes by 
taking the free axes of the observed points of error, or shot-marks. But if 
the axes in (43) are supposed to be rotated through the angle 0, and the 
coordinates to these new axes are x' and y' , so that we have 

x = x' cos — y' sin 6, y = x' sin 6 -\- y' cos 6, 
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then substituting these for x and y in (43), the exponent of e will contain 
not only x n and y' 2 , but the product x'y', so that it will not in general be 
possible to separate the second member into two factors, one containing the 
abscissa and the other the ordinate, as in (46). Thus the x and y deviations 
are not in general independent of each other unless the axes taken are free 
axes. But if we assume that the elementary errors are such as to make 
given deviations equally probable in all directions, every axis through the 
centre of gravity will be a free axis, and the deviations in horizontal and 
vertical directions will therefore be independent. That assumption is ordi- 
narily made in practice, as it simplifies the calculations. The deviations of 
the shot-marks are then measured from the centre of gravity directly, in- 
stead of from the axes. Denoting the square root of the mean of the squares 
of these deviations by edx, and taking dx = dy, we have 

s 2 = ef + e\, 
and since e x and e 2 are here equal, e 2 = 2ef = 2e|, and (44) gives 

(hydxf = (h 2 dyf = 1-W, .*. h x = h 2 = l-(edx). 
(Compare Didion's work already cited, and Sonnet's Dietionnaire des Math- 
e'matiques Applique es, article Probabilite du Tir.) The probability surface 
is what (43) becomes when dx = dy, h-^ = h 2 and x 2 -\-y 2 = r 2 , the elliptical 
sections being reduced to circles. It has been found to agree sufficiently 
well with the results of observation in target shooting. Probably it will not 
be worth while to compute the place of the free axes by (31) unless the ob- 
served points of error or shot-marks form a decidedly elongated group. If 
the probability of deviation is really the same in all directions, still it is not- 
likely that the actual distribution of the shots will be uniform, and a slight 
elongation of the group will not compel the inference that the elementary 
errors must be greater in that direction. 

Returning now to the strict construction of the surface (43), let us refer 
the position of the rectangle dxdy to polar coordinates by writing 

x — r cos cp , y = r sin <p, (47) 

so that we have 

__ h t h 2 dxdy -r*(h{ cos^+hlsm 2 <p) 

it 

as the probability that an error which occurs will fall on the area dxdy 

whose coordinates are r and <p . Suppose that the plane of XY is divided 

into rings by circles drawn from the origin as a centre, the breadth of each 

ring being j/(dxdy), and that it is also divided into sectors by radii drawn 

at angular intervals equal to d<p. The plane is thus partitioned into small 

elementary areas, each area at distance r from the origin being equal to rd<p 

X-\/(dxdy), and dxdy is contained in it rdf -f- ^/(dxdy) times, so that the 
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probability that an error which occurs will fall on this area, is 

zrdip-^-y '(dxdy). 
Hence the probability that it will fall on a given sector whose extreme 
radius is r will be 

P== ~dJ Z \ ] 7{d!cdy)) dr ' 
where dr is to be taken equal to \/{dxdy); so that substituting the value of 
z from (48), we have 

n J 
and by integration 

hlKdf_ __ J , -r^cos^+Afsin 2 ,*) \ ,. q . 

Making r = oowe get as the probability that the error will fall some- 
where on this sector extended to infinity, 

p = A 1 Vg___ (50) 

' 2w(AJ cos> + h\ sinV)' v ; 

If now it is required that p ~ %p, we have the relation 

-j- 2 (/i!cos 2 </>+/i,|sin 2 (&) 
e " =|, 

and restoring a; and 3/ from (47), 

*-«**+■**•> = x ... ^2 +A 2 y2 _ (i og 2)-*-(log e). 

Giving to Aj and A 2 their values as in (44), 

h\ = l-~2( ei dxf, hi = l+2(e 2 dyf, 

we get the equation of an ellipse, 

"* , I 2 _ -, 

1" r/i — a\ . n TUZ "X7T2 — x > 



[(log4)-4-(loge)](e 1< &)2 [(log4)--(loge)](e 2 rf2/) 2 
whose semi-axes are the square roots of the denominators in the first mem- 
ber. Their values are 

a = 1.17741(e 1 cfe), b = 1.17741(e 2 %), (51) 

so that they are proportional to e x dx and e 2 <%, the quadratic mean errors in 
the x and y directions. Hence for example, in the group of shot-marks al- 
ready considered, if we find the free axes and the quadratic mean deviations 
from them, and construct an ellipse with semi-axes a and b as in (51) coin- 
ciding with the free axes of X and Y respectively, this will be the ellipse of 
probable error, in the sense that it is an even chance whether a particular 
future shot will strike inside or outside the ellipse. And if two straight 
lines are drawn from its centre making any angle with each other, it is an 
even chance whether shots which fall between them will fall inside or out- 
side of the ellipse. 
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There is a correspondence between some of our results and those which 
have been reached by Arjdra and others in the theory of the error-ellipse 
(die Fehlerellipse), which is one of the more recent developm'ts of the method 
of least squares as applied to the location of a point by triangulations. The 
true position of the point should be at the common intersection of several 
straight lines, whose directions are given free from error, but whose positions 
laterally are subject to errors of observation, so that the lines do not in gen- 
eral have a common intersection. Each of them may be represented by an 
equation of the form 

ax -\- by -\- 1 = 0, 
where x and y are the coordinates of any point in the line, referred to rec- 
tangular axes taken in the vicinity of the desired point, a and — b are the 
sine and cosine of the acute angle which the line makes with the X axis, on 
the -fJTside, and I is the perpendicular from the origin upon the line. The 
angle and the perpendicular are essentially + or — according as thay lie 
above or below the Xaxis. The various intersecting lines furnish a number 
of observation equations of this form, which being solved by the method of 
least squares, give the desired coordinates x and y of the most probable 
point of common intersection. It is then shown that all other points whose 
probabilities are equal are ranged about the most probable point in ellipses 
which are similar and similarly situated and concentric. The directions of 
their major and minor axes are given by a formula which bears an external 
resemblance to (31). See for instance Jordan's Handbuch der Vermessungs- 
kunde, I., pp. 120 and 123. 

The dimensions of the ellipse within which it is an even chance that the 
true point will fall, have been ascertained by Helmert. (Ausghichungs- 
reohnung, p. 236.) He makes its semi-axes bear to the corresponding quad- 
ratic mean errors of the point, the same ratio 1.17741 which I have obtained 
in a different way in (51). Hence they bear to the corresponding probable 
errors the ratio 

1.1 7741 -=-.6745 =1.7456. 

His demonstration will be found in Sohlbmilch's Zeitsehrift fur Math, und 
Phys., 1868, pp. 75 to 79. The probability of a given deviation of the p't 
from its most probable place, is expressed by the same form of function 
which I have found as the equation (43) of the probability surface. 

It may be questioned, however, whether the theory of the error-ellipse 
affords proof of the universality of that law of facility, for all accidental er- 
rors in the situation of a point in a plane. It seems to me that when taken 
in this extended sense, the method is open to the same kind of theoretical 
objection which we noticed in connection with (46). An observed point of 
error, or intersection of two erroneous observed lines, is the result of the 
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concurrent happening of two independent errors, each of which, it is assum- 
ed by the theory, follows an exponential law of facility such as is expressed 
by the probability curve (1). The probability that any particular point in 
the plane is the true one, is regarded as proportional to the probabilty of the 
concurrent happening of all the independent errors represented by the per- 
pendicular distances from that point to all of the observed lines. The most 
probable point is that which renders the sum of the squares of the distances, 
each square being multiplied by the weight due to its observed line, a min- 
imum. (Helmert, in the Zeitschrift, pp. 76 and 89; also Jordan, II:, p. 223.) 
For the case of intersecting straight lines which lie in only two mutually 
perpendicular directions, the demonstration is evidently similar in principle 
to that which I remarked upon in connection with (46). My own proof of 
(43), however, is a direct proof, which does not take (1) into account at all, 
inasmuch as it does not assume that the x- and y- errors of a point happen 
independently of each other. But it includes the proof of (1) as a special 
case, for if the elementary points of error are all on one line which is made 
the axis of X, so that the given polynomial (2) is reduced to one dimension, 
all its coefficients being zero except those on the horizontal line through the 
middle of the square group, the last equation in (23) disappears, and of 
the two equations (30) only the first one remains, with a 2 = and y = 0. 
When the centre of gravity of the masses X is made the origin and centre of 
the linear group, by imagining terms with zero coefficients to be added if 
necessary, we get a 1 =0, and only the first of the two differential equations 
will exist in (35) and (39), the integral being 

The condition (41) is reduced to 



If*-', -, h,dx 

-j- I zdx — 1 , . ' . o = ', , 



and 2 becomes the ordinate to the probability curve, like y in (1). 

The occurrence of errors in two dimensions can be further illustrated by 
the drawing of balls from an urn, each ball being marked with two num- 
bers which may be positive or negative, the first number representing the x 
and the second the y coordinate of the point of error, reckoned from the 
true point, or place of error zero, as an origin. Each ball is replaced as 
soon as drawn, and they are all mixed so as to keep the chances the same. 

Suppose for example that the urn contains 20 balls marked as follows; 
one of them is marked 0, 0, four are marked 1, 0, five 1, 1, four 0, 1, one 
— 1, 1, three — 1, — 1, and two 0, — 1. Then if k drawings are made, the 
probability that the algebraic sum of the first or a-numbers on the balls 
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drawn will be t, and that at the same time the sum of their sec'd or y-num- 
bers will be v, is the coefficient of Ehf in the expansion of the polynomial 
(52) to the k power, this polynomial being a case under the general form (2). 

(52) 



m\ +1 +4f ) 



20 \+3e- 1 -^- 1 +2^- 

To find the centre of gravity and the free axes of the coefficients k, we 

may disregard for the present the divisor 20 common to them all, and write 

them in the rectangular form (53). Their centre of gravity rTHTg 

is at the point whose coordinates reckoned from the middle A -, \ , /KO v 
11 -> a , ■ . i a i a i 14 loo) 

term 1, are u' = \Au, v' = %Av, where Au and Av are the „,»„ 

intervals between the terms in (53), horizontally and verti- ' M l 
cally. Drawing horizontal and vertical axes of TJ and V through this cen- 
tre of gravity, the abscissas u of the terms in the three columns are 

— fz/w, — \Au, $ An, 

and the ordi nates v of the terms in the three rows are 

\Av, — \Av, —\Av. 

Denoting the terms in (53) by M, and applying (31), we find 
I\Muv) = 2^-AuAv, I(Mu 2 ) = ±?-{Au) 2 , 2{Mv 2 ) = $£-(Avf. 

^-«(s)*{«(5)'-«}- 

Taking for convenience Au = Av, so that (53) is a square, we have tan 2<p 
-— — $£-, wherefore ip = — 40°4'. This is the angle which one free axis 
makes with the horizontal, at the centre of gravity, and the second free axis 
is at right angles to the first. Designating these axes as X and Y, we get 
from (33) the relation between the new coordinates and the old, 

x = u cos — v sin 0, y = u sin + v cos 0, (54) 

where = — <p = 40°4'. Squaring x and y, multiplying them by M, and 
summing the results for all the terms M in the system, we have 

I(Mx 2 ) = 2'{Mu 2 ) cos 2 /? -f l\Mv 2 ) sin 2 /? — 22{Muv) sin cos 0, I ,,_. 
I{Mf) — Z{Mu 2 ) sin 2 /? + I{Mv 2 ) cos 2 /? + 2I(Muv) sin cos 0. / [i)0) 

We assign to the sums in the second members their values already found, 
and take for simplicity 

Ax = Ay = Au = Av, (56) 

and so get 

2\Mx 2 ) = 6.909(4c) 2 , ±Wy 2 ) = ISMS(Ay) 2 . 

Dividing these by 20, since X = M-t-20, we obtain the squares of the ra- 
dii of gyration of the masses k about the Fand Xaxes, 

r[ = .MM(Ax) 2 , r; = .9294(%) 2 . 
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Suppose that the number of drawings from the urn is k = 4. Then in the 
expanded polynomial the squares of the radii of gyration of the coeff's are 
hr\ = 1.382(Axf, h\ = 3.7l8(zfy) 2 , 

so that in (44) we have 

e\ = 1.382, e\ = 3.718, 
and (45) becomes 

z = .0702e- 862 ' 2 - 134 ^, \ ,. R . 

.-. log 3 = 2.8464— .1571 2i 2 — .05840/. J [ ] 

This is the equation of the limiting surface which represents approximately 
the contour of the coefficients in the expansion of (52) to the fourth power. 
Since the centre of gravity preserves in the expansion the same relative 
position which it had in the first power, its coordinates reckoned from the 
middle point will be four times as great, or 

u' = Au, v' = Av. 

The expansion contains 81 terms, wich occupy the centres of the small 
squares in the subjoined figure, and the centre of gravity falls at the next 
term diagonally from the middle one. The axes of U and V have their 
origin at this point, as do also those of X and Y. From the coordinates u and 
v for the several terms we compute the values of x and y by (54). A few 
nearest the origin being tabulated as in (59), the rest can be readily found 
from them by their differences. The values of z are now computed, and 
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being multiplied by 1000, they are as shown in the first half of the accompa- 
nying table, omitting decimals. Also by actual expansion of (52) to the 4th 
power we get the true coefficients, which being multiplied by 1000 like the 

(1000 z). 
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between the two sets of coefficients is evident, though it is not very close 
because h = 4 is not a large number. It will be seen that the largest tab- 
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ulated numbers form groups elongated in the direction of the free axis of Y. 
In accordance with (56), the numbers should properly be arranged in a 
square as in the figure. 

The position of the axes j 
of UV&nd of XY are there 
shown, and also the ellipse 
of probable error. The semi- 
axes of this ellipse as given | 
by (51) and (57), are 
a = l.'38dx, 
b ■= 2.27<%, 
where dx = dy is the side of I 
one of the small squares in 
the figure. 

The ellipse encloses rather I 
more than nine squares, al- 
lowing for fractions, but we I 
can see that it is not far from an even chance that the algebraic sum of the 
x- and y-numberson the balls drawn from the urn will be some one of the 
nine pairs. 

1, 2 2, 2 

1, 1 2, 1 

1, 2, 

By following the analogy of this ex- 
ample we can construct an equation such as (58) to represent the expansion 
of any given polynomial, and to show the distribution of the errors result- 
ing, in a given number of trials, from any set of elementary errors we please 
to assume. 

It should be noticed, however, that the coefficients of a polynomial repre- 
senting probabilities are all positive, whereas the general problem, to repre- 
sent the limiting form of an expansion when some of the coefficients are 
negative, may give rise to special cases not included in our present investi- 
gation. The limiting surface, if there is one, will in gen'l be the probab'ty 
surface (43), but not always. In (29) it may happen that /? 3 or /9 2 , or both 
of them, are zero, so that the differential equations (35) vanish, and our 
process fails. In such a case it is no longer sufficient to consider, as we did 
in connection with (26), that the element of the surface is a plane surface. 
Its curvature must be taken into account, and the limiting surface will 
most likely consist of an infinite series of undulations, analogous to those of 
the higher curves discussed in my former article. (See Analyst, Sept. '79) 



0,2 
0,1 
0,0 
The most probable result is 1, 1. 



